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In this article, we study the outage performance of an ambient multi-backscatter sys- 
tem supported by non-orthogonal multiple access (NOMA) downlink transmission, 
where legacy users receive information either directly from the base station (BS) or 
from multiple backscatter devices (BDs) as well. Specifically, we characterize the net- 
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work connections between the BS and network users into three modes of operation. In 
the first mode, the BS communicates with the users via the BDs by exploiting NOMA. 
In the second mode, the BS communicates directly with the users by employing selec- 
tion combining (SC) beamforming. In the third mode, the users receive signals from 
both the BS and BDs simultaneously, by using maximal ratio combining (MRC) to 
combine the BS and BDs links at the users. Consequently, we derive the exact outage 
expressions of each mode and utilize Monte Carlo simulations to validate the individ- 
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1. INTRODUCTION 

The rapid deployment of internet-of-things (loTs) in diverse systems such as energy systems (Smart 
Grid), healthcare systems, industrial systems, and transportation systems, has seen a tremendous push into re- 
searching ways of powering these billions of devices in a sustainable and environmentally friendly way [I], 
(2). Ambient backscatter communication (AmBC) has been proposed by the research community as a solution 
to address this unique challenge [3]. AmBCs can free up IoT devices from needing batteries as backscatter 
devices (BDs) harvest energy from traditional radio frequency (RF) transmitters such as television and radio 
transmitters to transmit information to nearby IoT devices [2]-{13]. Yang et al. in [£], the authors note that the 
harvested energy from ambient RF transmitters is adequate to provide power to high-rate battery-less sensors. 
Therefore, AmBC can be used in IoT environments where it’s not convenient to replace or maintain sensor bat- 
teries, e.g., in wide-body area networks (WBANs) [14], [15]. However, due to the broadcast nature of ambient 
RF transmitters, AmBC tend to suffer from direct-link interference from the RF transmitter [4]-[6]. Several 
methods such as the design of maximum-likelihood (ML) detectors, and interference cancellation techniques, 
have been proposed to eliminate the direct-link interference, see ial, i6] and references therein for more details. 

In this study, we propose utilizing non-orthogonal multiple access (NOMA) technology to address 
this problem, owing to its ability to simultaneously serve multiple wireless users over the same radio resource 
via linear superposition coding at the transmitter and successive interference cancellation (SIC) at the user’s 
receiver [7], [16)-[18]. Therefore, integrating NOMA with AmBC networks shows great promise for future 
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green IoT networks. Apart from NOMA, there have been several recent studies on deploying reconfigurable 
intelligent surfaces (RIS) to improve backscatter communication performance [19J—[21]. However, in this 
article, we focus on the contribution of NOMA to AmBCs performance. Our main contributions are: Firstly, 
we derive exact outage expressions for three practical AmBC link cases. In the first case, the base station (BS) 
communicates with the users via the BDs. In the second case, the BS communicates directly with the users. 
In the third case, the users receive signals from both the BS and BDs simultaneously; Secondly, from these 
expressions, we analyze the role of channel vectors, power allocation coefficients and backscatter reflection 
coefficient on outage probability. 

All expressions are confirmed using Monte Carlo simulations. The paper is organized is being as. In 
section 2, we describe the proposed system model and signal characteristics at the users. Thereafter, in section 
3, we derive exact outage probability (OP) expressions depending on the mode of operation. In section 4, we 
explain results, followed by the summary in section 5. 


2. RECEIVED SIGNALS AT TWO USERS 

In this paper, we study a downlink wireless system to allow the base station (BS) serving two multiple 
antennas users (user U; and user U2) under assistance of multiple backscatter devices (BDs), shown in Figure 
1. In some practical case, the second user U2 might connect with the BS directly. 


Figure 1. System model of multi-backscatter NOMA 


The criteria to choose best multi-backscatter is given by: 


x 2 
n* = ar max |h i 1 
me m On ( ) 


Peery 


BS broadcasts the information x?’ to multi-backscatter (B D»), n = 1,2,...N, which is given by: 


Ps = a,Psx1 + V/a2Psx2, (2) 


where xı and x2 are the messages with unit power transmitted to U and U2, respectively. Further, Ps is the 
total transmit power of the BS with power allocation parameters a, and ag, respectively, with aj > az and 
ag +a,=1. 

The maximal-ratio combining (MRC) is used with beamforming vector to achieve optimal transmis- 
sion scheme as in [22], with ||w;|| = 1: 


„iE {1, 2}, (3) 
where ||-|| denotes the Euclidean norm of a matrix. 


The BD the BS signal to U and U2 with its own message s, where E fls} = 1. Thus, U; and U2 


receive one type of signal: the backscatter link signal from the BD. The received signals at U; and U2 can thus 
be written as: 
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YU = Bho,n* hiw] rPSs EJ NU; (4) 
yu, = Bhon» ||h2wi]| £?Fs + nu, (5) 
where (@ is a complex reflection coefficient used to normalize, hı and ha are the 1 x M channel vector of 


BD — U; link and 1 x K channel vector of BD — U2 link, ny, and ny, are the additive white Gaussian noise 
(AWGN) with zero mean and variance matrix Tu, Im and Cu, I, respectively. Hence, the signal x, is detected 
successfully and x2 as interference. Then, the signal-interference to ratio at Uj is given by: 


pay B? \\hi ||" Rone 
Vii = 3 7 5 ; (6) 
pa2ß°||hil| |Ro,n*| +1 


where p = Ez, i € {1,2} is the transmit signal-to-noise radio (SNR). 
A 


Thus, the instantaneous SINR at U% to detect x; and the instantaneous SNR at U2 to detect x2 is given 
respectively as: 
2 2 2 
pa 8" || hall" [hon | 


2 2 2 
721 = r 72,2 = pa2b || ho|| [hon | : (7) 
pa23?||ha||"|hone|” +1 


3. OUTAGE PROBABILITY COMPUTATION OF DIFFERENT MODES 
3.1. Scheme 1: outage probability for backscatter-non-orthogonal multiple access system without direct 
link 
In this study, we assume independent Rayleigh-distributed random variables (RVs) for all the channels. 


Thus, RVs |ho,n«|"s 


n 


||h;||? have exponential distributions, with f n 2 (2) = DDAA A ROD! oT Pho, 
| O,n | n 


Qho 
ao ey 
N-1, py, x P-1 


final (x) = Eo and Fun, (x) =1- e Ohi 2 TEDO with T (x) = (a = 1)! is the gamma 


function. 
The OP of U; is: 


OPE, =Pr(n11 <1); (8) 


where 7; = 271 — 1 with Ry being the target rate at U; to detect x1. 
Theorem 1. The exact OP expression of U; is calculated as: 


M-1 N n—1 
N (a1 or 
OP, =1 -2 
Pū, D n ma 


m=0 n=1 


9 
6:9, \ 2 bın a 
x (32) Ki-m 2 5 
Qan Qn, Oha 
where 6, = Goa Note (9} is based on the condition of a; > 71a2 
Proof: 
OP, =1 — Pr | |\hil|? > — 
B?2|ho,n= | 
=1 F A t) dt 
51 J Fiat? (ga ) fron © (10) 
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_ co __ 91 _ nt 
where Final? (x) =1- Fin? (x) and ® (t) = J e P Aht Sho t”dt. 
Based on [24], (3.471.9). We can rewrite as: 


M-1 N n-1 
N (1) orn 
TN =]1—2 il 
OPi, >, > n JT (m+1) OF BPO, 


m=0 n=1 
an Lam J (11) 
144ho 1n 
x Ky_m | 2,4/ ——— |. 
(sos) i ( \ m) 
Next, the OP of U3 is given by: 
OPH, = Pr (422 < 12) = 1 — Pr (lhal? lho,n+|? > 62) (12) 


where y2 = 22 — 1 with Rə being the target rate at U2 to detect £2, 02 = = Similarly with solving OPS, 
can be achieved OPE, as: 


N n—lnk 
N (—1) ofn 
OPZ, =1-2 
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3.2. Scheme 2: outage probability of directed link by selection combining mode 
In the MRT scheme, to maximize the signal-to-interference noise (SINR) at U; and U2, BS uses a 


(13) 


beamforming transmit vector w, = ii , r € {SU,, SU>} to signal x?* before sending w, to users, hsuy, 
and Asu, are the 1 x M and 1 x K channel vector of BS to U; link and BS to U2 link, respectively. The BS 
transmits signals directly to two users as 


Yr = |hrwr|| £2? + np TE {SU1, SU2}. (14) 
We have SINR at Uj as: 


lsu: ll aip 
YSU, = 2 , (15) 
hsv, || a2p +1 
2 
and SINR and the successive interference cancellation (SIC) at Uz as ysu,, = paee and Uə as 
j SU a2p 
2 ; 
ysu, = ||hsu, || azp, respectively. 
Finally, the instantaneous SINRs at U; and U2 based on selection combining, can be written as: 
yS = Pr(max (ysu, Yii) < yi) iE {1,2} (16) 
The OP of U; is given by: 
= 2 x 
Fpa (0) = Pr (1al? < zir) 
== J Pima | (SEH) nont (t) dt (17) 
M-1 N N ee a 
—]— næ™(—1) 
=1 2 ( - ) TFI m Aq? (¢|z), 


where ô (t|x) = f e MSE Tho tm dt and A (x) = (a1 — raz) p8?. Using [24], (3.471.9) and after some 
0 


transform. F}, , (x) and F. 


ysu, (x) can be express respectively as: 


MAIN /N na™(—-1)"-1 
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(18) 
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and, 
M-1 Mo (ay zaz) egy 
ave 1 
F. =1 19 
sv, (2) mu T (m+ 1) gy, (a1 — zaz)” p” a 


Substituting and into (16), the close-form of OP of OPE is obtained as: 


M-1N nP (17-1 Nho l=m = 
= c n 
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M=1 gm, “pan (20) 
x |1 Š, TFD, 
Next, the OP of Uz is: 
OPi = Pr (3° < q2) = Pr (Ysu, < 72) Pr (q2,2 < %2). (21) 
Similarly, OPE can be achieved as: 
K-1 N l-k 
e EEC tn E] 


K-1 a 
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3.3. Scheme 3: outage probability of directed link by maximal ratio combining mode 
In this scheme, the backscattering link signals and direct link signals are combined by maximal ratio 
combining (MRC) at the two users. Hence, after using MRC the received SINR two users are given by: 


RO = ysu, + Yii i€ {1,2} (23) 


Then, the OP at U; can be calculated as: 


OPEET = Pr (ysu, < 71 — 11,1) 


V1 Yı— r (24) 
=f ft | Fa Wade. 
0 0 


It is very difficult to obtain the close-form solution of the integral of (24). Therefore, we use the 
approach in by making an I-step staircase approximation to the actual triangular integral region in (24) 
with sufficiently large I to obtain: 


= i+1 i 
aE (ru (2s) to (9) 
i=0 


(25) 
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Thus, the OP at U2 is given by: 
OPT = Pr (> MRC <72) 
02 0 
= F CEGA 
= J iseal (2)F li (3 =) (26) 
ax Ti- Tel 

= 1 2]; 

T(K) Rsv, 


Outage probability computation in multi-backscatter systems with multi-modes of operation (Dinh-Thuan Do) 


244 o ISSN: 2302-9285 


where 02 = oo Yı can be calculated as: 


2 eee, Sees 
1 K-1, %h — oK 0 
Yı = a e SU2 aan thee 5 (27) 
where y (x, y) is the lower incomplete Gamma function [24], (3.351.1). And YT. is written as: 


N K-i (k+1)/2(_ 1)” 1 
y y n 
vs z (A la 
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n A x 
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Putting t = 27/02 — 1 > 02 (t + 1)/2 = xz => 02/2dt = dz. But it is very difficult to obtain the close-form 
solution of T2, to do so, we use Gaussian-Chebyshev quadrature [26], (25.4.38), to get: 


To x 3 pe) ( r ) WAM Gt > ny IE EDET AECH? 


x 


sue 
N 


(k41)/2 G(F +1) 25K — K)Qe LDL 
nai Ko AA P J MPR Tees PO Oisy, PEHD) (29) 
nA (£) 
is, (2 ge). 


Substituting (26) and (28) into (25), OPE i is given by: 


K,02/ Qn sy, N K-1 L N 
OPEEE wy (Kite hsv) te) aS 2 n ) 
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where = cos (Alzin) and A (a) = (4 — SE. 


4. NUMERICAL RESULTS 

In this section, we set do = dı = dsu, = dsu, = 1, d2 = 2, Qn, = do", Qn, = di, Qn, = dy", 
Qsu, = dgy, and su, = dgy, and 7 is the path-loss exponent setting to be r = 2. In particular, main 
parameters can be seen in Table 1. In addition, the Gauss-Chebyshev parameter is selected as L = 100 to yield 
a close approximation. 


Table 1. System parameters used in the performance evaluation 


System Parameters Values 
Monte Carlo simulations repeated 10° iterations 
The power allocation coefficients {a1, a2} = {0.8, 0.2} 
The target rate at U1 Ry = 2 bps/Hz 
The target rate at U2 Rg = 2 bps/Hz 
The complex reflection coefficient B=0.5 


Figure 2 plots OP and base station transmit SNR relationship, with with varying N = kK = M = 1 
and N = K = M = 3. In (8), (12), (19), (20), (23), and (29) are used to plot the analytical lines. From Figure 
2, we notice that the network users experience different outage performances depending on the N = K = M 
values. The best performances are achieved by scheme 3, with the worst performance obtained by scheme 1. 
We can see the analytical curves fit with Monte Carlo simulations. 
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Figure 3 plots the OP relationship between outage and base station transmit SNR, with varying 
= 0.3 and 8 = 0.7. As in Figure 2, from Figure 3, we also observe that the users experience different 
outage performances depending on the complex reflection coefficient 3. The best performances are achieved 
by scheme 3, with the worst performance obtained by user 2 in scheme 1. 

Figure 4 plots the relationship between OP and yı = 7; (dB), varying 6 = 0.3, and 8 = 0.7. From 
Figure 4, we observe that some network users performance curves approach an outage probability ceiling. 
Generally, the outage performance reduces with increase in y; = 7; (dB). The best performances are achieved 
by scheme 3, with the worst performance obtained by scheme 1. 


—e— OP} - Ana. (6 

104 F| —7— OP iz" - Ana. 

=--O--- OP}, - Ana. 

=-= OPE, - Ana. 
r 


Outage Probability 
3 

Outage Probability 
3 


a 
+ 


-5 | |--A- OP EE - A 
ate OPHE - Ana, (8 = 0.7) 
---7--: OPH - Ana. (8 = 0.7) 
Sim. 


Figure 2. Outage probability versus p (dB) varying Figure 3. Outage probability versus p (dB) varying 
N=kK=M=1andN=kK=M=3 8 =0.3 and 8 = 0.7 


Ak] —o— OPE, - Ana. (8 =0.3) --++--- OPH, - Ana. (8 = 0.7) 
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—O— OPE - Ana. (G=0.3) ---A-- OPE - Ana. (8 = 0.7) 

—A— OPẸ - Ana. (6 =0.3) ---¥e-- OPE - Ana. (8 = 0.7) 
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yı =% (dB) 


Outage Probability 
3 
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Figure 4. Outage probability versus yı = yı (dB) varying 8 = 0.3 and 8 = 0.7 


Figure 5 plots the relationship between OP and power allocation coefficienct a2, varying N = K = 
M = 1 and N = K = M = 2 with p = 20 (dB) and Rı = Rə = 1. From Figure 5, we note that several 
OP curves approach an OP floor with increasing a2, while other OP curves converge at 0 when az = 1. Figure 
6 plots the relationship between OP and N = M = K, with p = 5 (dB). From Figure 5, we observe that 
maximal ratio combining delivers the best outage performance with the OP for backscatter-NOMA system 
performing the worst. In both Figures, user 2 in scheme | had the worst OP performance, due to the weak 
backscatter device signals that it recovers, therefore, it can easily experience outage in any condition. With 
user 1 in scheme 3 achieving the best performance due to the signals from the very near multi-backscattering 
devices and the further away base station direct link being combined by maximal ratio combining at its receiver. 
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Figure 5. Outage probability versus a varying N = K = M = 1 and N = K = M = 2 with p = 20 (dB) 
and Ry = Rə =1 


5. 


0 5 10 15 20 
N=M=K 


Figure 6. Outage probability versus N = M = K with p = 5 (dB) 


CONCLUSION 
In this paper, we provided the outage probability (OP) analysis of an ambient multi-backscatter system 


with multi-modes of operation. We derived exact solutions of OP for different users depending on the mode 
of operation. We observe from the simulation results, that the direct link between the base station and users 
affects the OP. In future work, we will consider ergodic capacity of the system. 
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